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\ Abstract 

The momentum density conjugate to a bosonic quantum field splits naturally into the sum of a 
■ classical component and a nonclassical component. It is shown that the field and the nonclassical 

\Q . component of the momentum density satisfy an exact uncertainty relation, i.e., an equality, which 

o : 

underlies the Heisenberg-type uncertainty relation for fields. This motivates a new approach to 

o . 

deriving and interpreting bosonic quantum fields, based on an exact uncertainty principle. In 

> \ 

particular, the postulate that an ensemble of classical fields is subject to nonclassical momentum 

D ! 

fluctuations, of a strength determined by the field uncertainty, leads from the classical to the 
.£h ' quantum field equations. Examples include scalar, electromagnetic and gravitational fields. For 

X ■ 

the latter case the exact uncertainty principle specifies a unique (non-Laplacian) operator ordering 

a ; 

for the Wheeler-deWitt equation. 
PACS numbers: ll.10.Ef, 03.70.+k 
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I. INTRODUCTION 



The Heisenberg uncertainty relation 



AxAp >h/2 



(1) 



for the rms position and momentum uncertainties of a quantum particle is a well known 
feature of quantum mechanics. It has recently been shown, however, that this relation is 
a consequence of a more fundamental connection between the statistics of complementary 
quantum observables. 

In particular, the distinguishing "nonclassical" property of complementary observables 
is that they cannot be simultaneously measured to an arbitrary accuracy. It is therefore 
natural to consider the decomposition of one such observable, the momentum say, into the 
sum of a "classical" and a "nonclassical" component: 



where the classical component, p cl , is defined as that observable closest to p (in a statistical 
sense) which is simultaneously measurable with the complementary observable x [1]. 

It turns out that such decompositions do indeed, in a number of ways, neatly separate clas- 
sical and nonclassical properties of quantum observables. For example, for one-dimensional 
particles the nonclassical component of the momentum satisfies the uncertainty relation 



for all pure states, where Sx denotes a measure of position uncertainty from classical statistics 
called the Fisher length [1, 2]. This exact uncertainty relation is far stronger than (and 
implies) the corresponding Heisenberg uncertainty relation in Eq. (1). 

The surprising fact that quantum systems satisfy exact uncertainty relations has recently 
provided the basis for deriving much of the quantum formalism, from an exact uncertainty 
principle. In particular, the assumption that a classical ensemble is subjected to nonclassical 
momentum fluctuations, of a strength inversely proportional to uncertainty in position, has 
been shown to lead directly from the classical equations of motion to the Schrodinger equa- 
tion [3]. A brief overview of exact uncertainty properties of quantum particles is provided 
by the conference paper in Ref. [4]. 



p = p + p' 



JIC 



5xAp nc = 



h/2 



(2) 
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The aim of this paper is to show that the abovementioned results can be generalised 
to bosonic quantum fields. The paper is therefore divided into two separate and logically 
distinct sections: exact uncertainty relations in Sec. II, and an exact uncertainty principle 
in Sec. III. Each of these sections may be read without reference to the other, according to 
the reader's interest. 

In Sec. II it is shown that the statistics of any bosonic field / and its conjugate momentum 
field g are connected by an exact uncertainty relation analogous to Eq. (2). This uncertainty 
relation provides a precise equality between the statistical covariance of the non-classical 
component of g and the Fisher covariance of /, and underlies a Heisenberg-type inequality 
relating the statistical covariances of / and g. 

In Sec. Ill it is shown that the quantum field equations follow from an exact uncertainty 
principle, for the case of bosonic fields with Hamiltonians quadratic in the field momenta 
(eg, scalar, electromagnetic and gravitational fields). This "exact uncertainty" approach is 
extremely minimalist in nature: unlike canonical quantisation, it does not use nor make any 
assumptions about the existence of operators, Hilbert spaces, complex amplitudes, inner 
products, linearity, superposition, or the like. The sole "nonclassical" element needed is 
the addition of fluctuations to the momentum density of members of a classical ensemble of 
fields, with the fluctuation statistics assumed to be determined by the ensemble statistics. 
The exact uncertainty approach is thus conceptually very simple, being based on the core 
notion of statistical uncertainty (intrinsic to any interpretation of quantum theory). 

As a bonus, the exact uncertainty approach further implies a unique operator ordering 
for the quantum field equations - something which the canonical quantisation procedure 
is unable to do. This is relevant in particular to the Wheeler-deWitt equation (where the 
associated unique ordering is, moreover, consistent with Vilenkin's "tunneling" boundary 
condition for inflationary cosmology [5]). 

It is remarkable that the basic underlying concept - the addition of "nonclassical" mo- 
mentum fluctuations to a classical ensemble - carries through from quantum particles to 
quantum fields, without creating conceptual difficulties. This is a strength of the exact 
uncertainty approach, not mirrored in other approaches that rely on connecting the equa- 
tions of motion of classical and quantum ensembles. For example, the so-called "causal" 
interpretation of Bohm and co-workers is explicity non-local [6] (and hence non-causal for 
relativistic fields), while one cannot simultaneously define both the electric and magnetic 
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fields in generalisations of Nelson's stochastic approach to electromagnetic fields [7]. 

Necessary elements from classical field theory and functional analysis are briefly sum- 
marised in the Appendices. 

II. EXACT UNCERTAINTY RELATIONS FOR FIELDS 

A Heisenberg-type uncertainty relation may be given for scalar fields, in the form of an 
inequality for the covariance functions of the field and its conjugate momentum density. 
However, a stronger result follows via a natural decomposition of the momentum density 
into classical and nonclassical components, and the use of the Fisher covariance function 
from classical statistical estimation theory. In particular, an exact uncertainty relation for 
scalar fields is obtained in Sec. II. D below, which underlies and implies the corresponding 
Heisenberg-type uncertainty relation, and which may be generalised to arbitrary bosonic 
fields. 

This derivation of exact uncertainty relations, i.e., precise quantitative connections be- 
tween the statistics of conjugate fields, from the quantum formalism, is to some extent 
conceptually reversed in Sec. Ill, where the quantisation of a large class of classical fields 
is achieved via the use of an "exact uncertainty" principle. The reader primarily interested 
in this new approach to obtaining the quantum field equations may wish to skip directly to 
Sec. III. 

A. Heisenberg-type relations 

We begin with the simplest case: a real scalar field / with conjugate momentum density 
g. Spatial coordinates will be denoted by x (irrespective of dimension), and the values of / 
and g at position x by f x and g x respectively. 

The underlying origin of the field is not at issue here - it could be a relativistic field 
satisfying the Klein-Gordon equation [8] , or a nonrelativistic field describing vibrations of a 
stretched string [9]. What is important for the purposes of this Section is that the corre- 
sponding quantised fields / and g satisfy the equal-time commutation relations [8] 

[fx, fx'] = [9x, 9x'\ = 0, [f x , g x '\ = iK8(x - x'). (3) 
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Thus, in the Schrodinger representation, where the quantum state of the field at a given 
time is formally described by a complex amplitude functional ^[/], the action of / and g is 
given by [10] 

= ^ = - 77 , (4) 

where 5/5f denotes the usual functional derivative (see also Appendix A). 

The non-vanishing commutator in Eq. (3) reflects the complementary nature of / and g: 
in general one must choose to accurately measure either the field or its momentum density. 
This is analogous to the complementary nature of the position and momentum observables 
X and P of a system of quantum particles, with non- vanishing commutator [Xj, P k ] = ih5j k . 
In the latter case one has the corresponding Heisenberg uncertainty relation [11] 

Cov(X) Cov(P) > (h 2 /4)I (5) 

for the covariance matrices of X and P, where [Cov(A)]j fc := (AjA k ) — (Aj)(A k ), and / 
denotes the identity matrix. Thus, since the covariance matrix vanishes for the dispersion- 
free case, the complementary observables X and P cannot be simultaneously specified. 
The covariance function of a given field h is defined by 

[Cov(h)] xx , := {h x h x i) - (h x )(h x ,), (6) 

where ( ) denotes an ensemble average. Thus the diagonal elements of the covariance function 
(i.e., x = x') represent the variance of the field at each point, and the off-diagonal elements 
correspond to the degree of linear correlation between field values at different points. 

The covariance functions of the conjugate fields / and g satisfy the Heisenberg-type 
inequality 

Cov(/) * Cov(#) y (h 2 /A)t. (7) 

analogous to Eq. (5), where we use a matrix-type notation with multiplication A* B defined 
by 

i A * B \xx> '■= J dx" A xx „B x n x ^ 
multiplicative identity 1 defined by 

t xx , := 5(x - x'), 
and the ordering A y B holds if and only if 

rf(A - B)h := J dxdx' (A - B) xx ,h* x h x , > 
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for all h. Eq. (7) is clearly a continuous analogue of the Heisenberg inequality in (5), and 
indeed may be "proved" from the latter inequality by a discretization argument (where the 
spatial degree of freedom is represented by discrete cells of volume e at positions Xj, the 
identifications Xj = f x , Pj = eg Xj , eJ2j = f dx are made, and the limit e — > is taken). A 
more rigorous proof is given in Sec. II. D below. 

Eq. (7) represents a field-theoretic generalisation of the 1-dimensional Heisenberg un- 
certainty relation in Eq. (1). Moreover, it has the same physical significance: since the 
covariance function of field h vanishes for the dispersion- free case, the conjugate fields / and 
g cannot simultaneously be specified. It will be shown in Sec. II. D that this inequality can 
be replaced by an exact uncertainty relation analogous to Eq. (2). However, two concepts 
must first be introduced: decomposition of the quantum momentum density into "classical" 
and "nonclassical" components, and the Fisher covariance function. 

B. Decomposition of the momentum density 

Since, in contrast to classical fields, the quantum momentum density cannot be measured 
simultaneously with the field, we define the classical component of the momentum density 
to be that observable which is closest to g, statistically speaking, under the constraint of 
being comeasurable with /. Thus measurement of the classical component provides the best 
possible estimate of the momentum density compatible with a measurement of the field. 

More formally, let g cl denote the classical component of the momentum density g. The 
comeasurability of / and g cl then implies that the action of the latter in the representation 
of Eq. (4) has the form 



for some real functional g d [f]. Further, for the classical component to be as close as possible 
to the quantum momentum density, the average deviation of g cl and g must be as small as 
possible, i.e., 



at each position x. 

Conditions (8) and (9) determine the classical component of the momentum density 



g d *=g d [f]V 



(8) 




(9) 
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uniquely, for each quantum state one finds, as shown further below, that 

dr „. h \l 5® 1 5**1 

* [/] = 2^57-^TrJ- (10) 

Note that the classical component depends on ^ - to make the best possible estimate one 
needs some knowledge about the state of the field [12]. Thus one should, strictly speaking, 
use a notation such as g d '^ for the classical component, particularly if one wishes to calculate 
expectation values such as (<7 d '*)$ for states $ other than \I>. However, since quantities 
involving the classical component will in fact only be evaluated for the corresponding state 
\& in what follows, such explicit notational dependence on the state may be conveniently 
dispensed with. 

Note also (for each state ty) that, since the operator g d commutes with the field operator 
/, the former operator is measurable to the same degree as the latter. In particular, the 
experimentalist's in-principle procedure for measuring g cl on state \1/ is to (i) prepare the 
quantum state (ii) measure the field /; and (iii) for measurement result / for / calculate 
the corresponding measurement result g cl [f] for g cl . 

It follows that the momentum density of a quantum field admits a natural (state- 
dependent) decomposition, 

9 = 9 d + 9 nc , (11) 
into the sum of a classical component g cl , corresponding to that part of g which is comeasur- 
able with /, and a nonclassical component g nc , corresponding to an intrinsically "quantum" 
remainder. Further, these components satisfy the relations 

<$*> = <$£>, <£T>=o, (12) 

(A^) 2 = (Agff + (A^ c ) 2 , (13) 
and hence the classical and nonclassical components also represent an "average" part and a 
"fluctuation" part respectively, the uncertainties of which add in quadrature. 

To derive the explicit form in Eq. (10) for the classical component, consider any other 
field operator h comeasurable with /, so that h^/ = h[f]^f in analogy to Eq. (8). Then from 
Eqs. (4) and (10), and Eq. (A3) of Appendix A (assuming (h) is finite), one has 

h 



h f 

(g x h x + h x g x ) = - J Df * 



S(h x V) 



8 fx 8 fx 

(h/i) J Df [**(6*/6f x ) - ^(5^*/5f x )]h a 

2(9fhx), 



where (A) is defined as the functional integral / Df^*A^ (see Appendix A). Hence 

<(& - kf) = ((^) 2 ) + ((k) 2 } - 2(gfK) 

Since the last term on the right is nonnegative, the lefthand side is minimised for the choice 
h = g d , as required by condition (9). 

Note that Eq. (12) follows immediately from Eqs. (10) and (11), and Eq. (13) then follows 
by substituting h = g cl into the expression immediately above. Eq. (13) in fact corresponds 
to the diagonal elements of the more general covariance relation 

Cov(g) = Cov(g d ) + Cov(g nc ), (14) 

which may be derived by explicit calculation of Cov(g — g cl ). The classical and nonclassical 
field components are therefore linearly uncorrelated. 

C. Fisher covariance 

Many different classical measures of statistical uncertainty, such as variance and entropy, 
have been used in writing Heisenberg-type inequalities for quantum particles [13]. Similarly, 
there is a wide degree of freedom in choosing measures of uncertainty for quantum fields - the 
covariance function Cov(f) appearing in Eq. (7) only represents one particular choice. An 
alternative choice is provided by classical statistical estimation theory: the Fisher covariance 
function. 

The Fisher covariance function Covp(/) is, like the standard covariance function Cov(/), 
a positive and symmetric two-point function which vanishes for dispersion-free ensembles, 
and is hence a suitable candidate for describing field uncertainty over an ensemble. Its main 
role in classical statistics is as the "minimum" covariance function associated with estimates 
of translations of the field. 

To motivate Covp(/), consider a classical ensemble of fields described by some probability 
functional P[f], with ensemble average (/) = / Df Pf = 0. If each member of the ensemble 
is translated by some amount /, a new ensemble is obtained, described by P[f — f] and 
having field average (/) = /. The two ensembles thus differ only by the displacement field 
/, and hence distinguishing between two such ensembles is equivalent to estimating this 
relative displacement. 
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How easily can such ensembles be distinguished? The basic answer is that any unbiased 
estimate / of the displacement field (i.e., with (/) = /) has a covariance function greater 
than or equal to the Fisher covariance: 

Cov(/) h Cov F (/). (15) 

The Fisher covariance therefore characterises the minimum uncertainty of such estimates. 
For example, the quantity 

(A F f x f : = [Cov F (/)L 

represents a lower bound for the variance, (Af x ) 2 , of any unbiased estimate of the displace- 
ment field at position x. Roughly speaking, two ensembles differing by displacement field 
/ stand a decent chance of being distinguished at position x if the Rayleigh-type criterion 
|/ x | > A F f x is satisfied. 

To explicitly define the Fisher covariance function, let F(f) denote the "Fisher infor- 
mation matrix" of an ensemble described by probability density functional P[f] [11, 14], 
with 

Then Covf(/) is defined as the corresponding matrix inverse, i.e., 

Cov F (/) * F(f) = F(f) * Cov F (/) = 1. (17) 

Fisher information was introduced in 1925 for one-dimensional statistical variables [15], and 
Eq. (15) is the continuous matrix version of the Cramer- Rao inequality [14], well known in 
classical statistics. A particular case corresponds to the choice / = / in Eq. (15), yielding 

Cov(/) h Cov F (/). (18) 

The Fisher covariance is therefore also a lower bound for the standard covariance. As this 
result is needed for establishing the connection between Heisenberg and exact uncertainty 
relations in Sec. II. D, it is useful to give a direct proof here. 

First, let A[f] denote an arbitrary functional of / with (A) < oo, and define 
/[5\nP 5A\f5\nP 5A\\ 

5 A 5 A \ r _ fSP 5A 5P 5A\ 



l dA dA \ r IdfdA 
[ U)Lx ' + \SKsL'/J WM + 6f x ,6f x 



rrn/ „ N1 ISA 5 A 5 2 A 



8 fx S fx' Sf x 5f x 



Restricting A to have the quadratic form 



dxdx' K xx ,(f x - h x )(f xl - h x >) 



(20) 



for some symmetric matrix function K and field h, and noting that M >z from Eq. (19), 
then gives 



where we have chosen h = (f) at the end of the calculation. Finally, Eq. (18) follows by 
choosing K to be the matrix inverse of Cov(f). Note that equality holds in Eq. (18) if and 
only if M = in Eq. (19), i.e., if and only if S(]nP)/Sf = SA/5f. Thus this Cramer-Rao 
inequality is saturated if and only if P ~ exp(A), i.e., noting Eq. (20), if and only if the 
ensemble is Gaussian. 

It is important to emphasise here that the Fisher covariance function is a purely classical 
quantity, requiring no reference to quantum theory whatsoever for its motivation or defini- 
tion. However, it plays an important role as a measure of uncertainty in an exact uncertainty 
relation for quantum fields, to be derived in Sec. II. D. 

D. An exact uncertainty relation 

It has been shown that the momentum density of a bosonic quantum field admits a 
natural decomposition into classical and nonclassical components. Moreover, Eqs. (3), (8) 
and (11) imply that 



Hence it is the nonclassical component which is responsible for the non-vanishing commuta- 
tor in Eq. (3). However, this nonclassical component does not merely satisfy a Heisenberg- 
type inequality analogous to Eq. (7): it satisfies the exact uncertainty relation 



Thus, the field uncertainty precisely determines the nonclassical momentum uncertainty, and 
vice versa. 

The exact uncertainty relation in Eq. (21) is the main result of this Section. It is 
remarkable that natural measures of uncertainty can be chosen which precisely quantify the 
uncertainty principle for bosonic fields. Note that the quantities appearing in Eq. (21) above 



F(f) y 2K - K * Cov(f) * K, 



[Lg c J} = o, {Lg:?}=ih5(x-x>). 



Cov F (/) * Cov(g nc ) = (h 2 /A)l. 



(21) 
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are operational in nature: Covf(/) may be calculated from the statistics of /, and Cov(g nc ) 
may be calculated from the statistics of / and g, via Eq. (14). Thus the exact uncertainty 
relation directly relates the statistics of the conjugate fields / and g. 

To demonstrate Eq. (21), note that for an ensemble of quantum fields described by 
probability functional P[f] = \^[f}\ 2 , the Fisher information follows from Eq. (16) as 



+ 2 Df 



T 1^ 

OJx 

* T7~ c - C - 

OJx 

+ c.c. 



T-l M 

* 71^ + C.C. 
OJx' 



8 fx 8 f x 

= -(4/R 2 )<^> + (4/R 2 )^>, 
where the last line follows via Eqs. (4) and (10). Hence, using Eqs. (12) and (14), 

CovGT) = (h 2 /A)F(f), 



(22) 



i.e., i/ie covariance of the nonclassical momentum density is proportional to the Fisher in- 
formation of the field. Eqs. (21) and (22) are equivalent. 

The exact uncertainty relation in Eq. (21), being a strict equality for all states of the field, 
is clearly far stronger in character than the corresponding Heisenberg-type relation in Eq. (7). 
Moreover, it immediately implies the latter, since Cov(g) >z Cov(g nc ) and Cov(/) >z Cov F (/) 
from Eqs. (14) and (18) respectively. 

Various applications of such exact uncertainty relations may be made, analogous to those 
made in Refs. [1,2] for quantum particles, but will not be pursued here. However, the notion 
that the momentum density of a bosonic field decomposes into a classical and a nonclassical 
part, with the uncertainty of the latter determined by the uncertainty of the field, underlies 
the formulation of a new approach to quantum field theory in Sec. Ill below. 



E. Multicomponent and complex bosonic fields 

It is of interest to briefly outline the form of exact uncertainty relations for more gen- 
eral bosonic fields. Thus, let f a denote a multicomponent bosonic field with conjugate 
momentum density g a . For example, f a may be a complex Klein-Gordon field <j), or the 
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electromagnetic field A^. This notation also covers the case where the index a labels several 
interacting bosonic fields. 

The non-vanishing equal-time commutation relation is now 

[/« =ih5 ab 5(x-x'), (23) 

and the momentum density operator decomposes into classical and nonclassical components 
as before, with 

g a = g cl > a + (f c ' a ? (24) 



(25) 



5f- v ' 6f° 

Clearly the components of g d commute with each other; it may also be checked that the 
components of g nc commute (which is necessary for the covariance function Cov(g nc ) to be 
well defined). 

The generalised covariance and Fisher information matrices are defined by 



[Cov(/o]£ :=((Kyh b x ,)-((Ky)(h b x! ), (26) 

in analogy to Eqs. (6) and (16) respectively, and similarly are purely classical statistical 
quantities. 

Finally, the corresponding exact uncertainty relation has precisely the form of Eq. (21), 
providing that the definitions of the multiplicative identity and matrix multiplication are 
replaced by the natural generalisations 



lA*BC:=Y,fdx"A:i,,BX>- 



respectively. 
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III. QUANTISATION OF FIELDS FROM AN EXACT UNCERTAINTY PRINCI- 
PLE 



It has been shown that the momentum density of a bosonic quantum field splits naturally 
and neatly into the sum of classical and nonclassical parts, as per Eq. (11). Moreover, the 
uncertainty of the nonclassical part is precisely linked to the uncertainty of the field itself, 
as per Eq. (21). These results not only provide a mathematical link between the statistics 
of conjugate quantum fields: they suggest the existence of a new bridge connecting classical 
field theory to quantum field theory - an exact uncertainty principle. 

Indeed, as demonstrated in this section, the assumption that an ensemble of classical 
fields is subjected to "nonclassical" momentum fluctuations, of a magnitude determined by 
the field uncertainty, leads from the classical equations of motion to the equation of motion of 
a bosonic quantum field. This result is restricted to the case of field Hamiltonians quadratic 
in the momentum density, but nevertheless covers most classical fields of physical interest, 
including the Klein-Gordon, electromagnetic and gravitational fields. 

A. Classical ensembles 

We first consider a real multicomponent classical field / = (f a ) with conjugate momentum 
density g = (g a ), described by some Hamiltonian functional H[f,g,t]. Spatial coordinates 
are denoted by x (irrespective of dimension), and the values of field components f a and g b 
at position x are denoted by /° and g b x respectively. 

The equations of motion for an ensemble of such fields are given by the Hamilton- Jacobi 
equation 

(99 

— + H[f,6S/6f,t} = 0, (28) 

and the continuity equation 




for the momentum potential S and the probability density P, as is reviewed in Appendix 
B. These equations specify the motion of the ensemble completely, where the momentum 
density associated with / is given by g — 5S/5f. 
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Defining the ensemble Hamiltonian H as the functional integral 

H[P, S, t] := (H) = jDf PH[f, 5S/5f, t], (30) 

the equations of motion for the dynamical variables P and S may be written in the Hamil- 
tonian form 

dP_ _ 6H_ dS_ 6H_ 

~dt ~ US' ~dt~ ~JP' ( ' 

The variational derivative of a functional integral such as H, with respect to a functional 
such as P or S, is discussed in Appendix A, and the equivalence of Eqs. (28) and (29) to 
Eqs. (31) follows directly from Eq. (A6). Hence, in analogy to Eqs. (Bl) of Appendix 
B, P and S may be regarded as canonically conjugate junctionals, which evolve under the 
ensemble Hamiltonian H [16]. Note from Eq. (30) that H typically corresponds to the mean 
energy of the ensemble. 

In what follows, we will only consider classical fields with Hamiltonian functionals 
quadratic in the momentum field density, i.e., of the form 

H[f,g,t] =J2j dxK?[f]g a x g b x + V[f}. (32) 

a,b 

Here K x b [f] = K ba [f] is a kinetic factor coupling components of the momentum density, and 
V[f] is some potential energy functional. The corresponding ensemble Hamiltonian is given 
by Eq. (30). Note that cross terms of the form g%g h x , with x ^ x' are not permitted in local 
field theories, and hence are not considered here. 



B. Momentum fluctuations =>• quantum ensembles 

The ensemble Hamiltonian H[P,S,t] in Eq. (30) is our classical starting point for de- 
scribing an ensemble of fields. This starting point must be modified in some way if one 
is to obtain new equations of motion, to be identifed as describing a quantum ensemble of 
fields. For example, the standard approach (canonical quantisation) assumes the existence 
of a complex field functional ty[f, t] living in some Hilbert space, and one simply writes 
down the Schrodinger equation ihd^/dt = H[f,—ih(5/5f),t]^ [10]. While this approach 
can generally be made to work successfully, it provides no explanation for the appearance of 
the nonclassical objects \& and h, a Hilbert space, and a linear operator equation. Moreover, 
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for cases in which K® h in Eq. (32) is dependent on / (eg, gravitational fields), the canonical 
approach is ambiguous as to the ordering of the functional derivative operator. 

In contrast, we take a very different approach, where the only a priori nonclassical feature 
is an assumption that the classical ensemble Hamiltonian H must be modified to take into 
account the existence of nonclassical fluctuations of the momentum density, with the magni- 
tude of the fluctuations determined by the uncertainty in the field. This "exact uncertainty" 
approach is motivated by the fact that bosonic fields satisfy exact uncertainty relations such 
as Eq. (21) of Sec. II. It leads to equations of motion equivalent to the Schrodinger equation 
above, with the added advantage of a unique operator ordering. 

Suppose then that SS/ Sf is in fact an average momentum density associated with field 
/, in the sense that the true momentum density is given by 



where N is a fluctuation field that vanishes on the average for any given field /. Thus 
the physical meaning of S changes to being an average momentum potential. No specific 
underlying model for N is assumed or necessary: in the approach to be followed, one may in 
fact interpret the "source" of the fluctuations as the field uncertainty itself. Thus the nature 
of the fluctuation field is not important - its main effect is to remove any deterministic 
connection between / and g. 

Since the momentum fluctuations may conceivably depend on the field /, the average 
over such fluctuations for a given quantity A[f, N] will be denoted by A[f], and the average 
over fluctuations and the field by (A). Thus N = by assumption, and in general (A) = 
J Df P[f]A[f]. Assuming a quadratic dependence on momentum density as per Eq. (32), it 
follows that when the fluctuations are significant the classical ensemble Hamiltonian H = 
{H) in Eq. (30) should be replaced by 



a,b 

Thus the momentum fluctuations lead to an additional nonclassical term in the ensemble 
Hamiltonian, dependent on the covariance matrix Cov x (iV) of the fluctuations at position 



g = ss/Sf + N , 



(33) 



H' = (H[f,5S/5f + N,t] ) 




(34) 



15 



x, where 

[Cov x (N)] ab :=N*Nl. (35) 

The additional term in Eq. (34) is uniquely specified, up to a multiplicative constant, by 
the following four assumptions: 

(1) Causality: H' is an ensemble Hamiltonian for the canonically conjugate functionals P 
and S, which yields causal equations of motion. Thus no higher than first-order functional 
derivatives can appear in the additional term, implying that 

Cov x (N) = a(P, 5P/Sf x , S, 5S/Sf x , f x , t) 

for some symmetric matrix function a. 

(2) Independence: If the system comprises two independent non-interacting subsystems 
1 and 2, with factorisable probability density functional P[f (1 \ f {2) ] = Pi[/ (1) ]P 2 [/ (2) ], then 
the dependence of the subsystem fluctuations on P only enters via the corresponding prob- 
ability densities P\ and P<i respectively. Thus 

Cov x (NW)\ pip2 = Cov x (NW)\ pi , Cov^ 2 ))^ = Cov x (N^)\ p2 

for such a system. Note that this assumption implies that ensemble Hamiltonians are addi- 
tive for independent non-interacting ensembles (as are the corresponding actions [16]). 

(3) Invariance: The additional term transforms correctly under linear canonical trans- 
formations of the field components. Thus, noting that / — > A -1 /, g — > A T g is a canonical 
transformation for any invertible matrix A (with transpose A T and coefficients A^), which 
preserves the quadratic form of H in Eq. (32) and leaves the momentum potential S invari- 
ant (since 8/8 f — > A T 5/5f), one has from Eq. (33) that — > A T A^, and hence that 

Cov x (N) -> A T Cov x (N)A for / -> A" 1 /. 

Note that for single-component fields this reduces to a scaling relation for the variance of 
the fluctuations at each point x. 

(4) Exact uncertainty: The uncertainty of the momentum fluctuations at any given po- 
sition and time, as characterised by Cov x (N), is specified by the field uncertainty at that 
position and time. Thus, since the field uncertainty is completely determined by the prob- 
ability density functional P, it follows that Cov x .(A^) cannot depend on S, nor explicitly on 
t. 
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It is seen that the first three assumptions are essentially classical in nature, while the 
fourth assumption is not: it postulates an exact connection between the nonclassical mo- 
mentum uncertainty and the field uncertainty. Remarkably, these assumptions lead directly 
to equations of motion of a bosonic quantum field, as shown by the following Theorem and 
Corollary (proofs are given in Appendix C). 

Theorem: The above assumptions of causality, independence, invariance, and exact 
uncertainty imply that 

Apf = C(5P/5r x )(5P/5f b x )/P 2 , (36) 

where C is a positive universal constant. 

The theorem thus yields a unique form for the additional term in Eq. (34), up to a 
multiplicative constant C . The classical equations of motion for the ensemble are recovered 
in the limit of small fluctuations, i.e., in the limit C — > 0. Note that one cannot make the 
identification N x ~ (5P/5f£)/P from Eq. (36), as this is inconsistent with the fundamental 
property N% = 0. 

The main result of this Section is the following corollary (proved in Appendix C): 
Corollary: The equations of motion corresponding to the ensemble Hamiltonian H' can 
be expressed as the single complex equation 

l h^ = H[f,-th5/5f,t]^ = -h 2 / dx^K?\f\-^J* + V\f\V, (37) 

where one defines 

h := 2VC, := \fPe iS/n . (38) 

Eq. (37) may be recognised as the Schrodinger equation for a quantum bosonic field, and 
hence the goal of deriving this equation, via an exact uncertainty principle for nonclassical 
momentum fluctuations acting on a classical ensemble, has been achieved. The ensemble 
of fields corresponding to ensemble Hamiltonian H' will therefore be called the quantum 
ensemble corresponding to H. It is remarkable that the four assumptions of causality, 
independence, invariance and exact uncertainty lead to a linear operator equation. 

Note that the exact uncertainty approach specifies a unique operator ordering, 
(8/8f x )K x lb (5/5fl), for the functional derivative operators in Eq. (37). Thus there is no 
ambiguity in the ordering for cases where K x b depends on the field /, in contrast to tradi- 
tional approaches (eg, the Wheeler-deWitt equation, discussed in Sec. III.D below). The 
above results generalise straightforwardly to complex classical fields. 
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Finally, it may be remarked that the equations of motion of a classical ensemble may be 
subject to some imposed constraint(s) on /, P, S and H. For example, each member of 
an ensemble of electromagnetic fields may have the Lorentz gauge imposed (see Sec III.C 
below). As a guiding principle, we will require that the corresponding quantum ensemble is 
subject to the same constraint (s) on /, P, S and H' . This will ensure a meaningful classical- 
quantum correspondence for the results of field measurements. However, consistency of the 
quantum equations of motion with a given set of constraints is not guaranteed by the above 
Theorem and Corollary, and so must be checked independently for each case. 

C. Example: Electromagnetic field 

The electromagnetic field is described, up to gauge invariance, by a 4-component field 
A^. In the Lorentz gauge all physical fields satisfy d^A^ = 0, and the classical equations of 
motion in vacuum are given by d v d v A IJj = 0. These follow, for example, from the Hamiltonian 
[10] 

H GB [A, tt] = (1/2) J dx V , u (tt'V - VA» ■ VA») , (39) 

where r]^ denotes the Minkowski tensor, 7r M denotes the conjugate momentum density, 
and V denotes the spatial derivative. Here Hqb corresponds to the gauge-breaking La- 
grangian L = —(1/2) / dx A^ v A^ v , and is seen to have the quadratic form of Eq. (32) (with 
K£ v = rj^/2). Hence the exact uncertainty uncertainty approach immediately implies that 
a quantum ensemble of electromagnetic fields obeys the Schrodinger equation 

ih(dV/dt) = H GB [A, -ih(5/5A)]V, (40) 

in agreement with the Gupta-Bleuler formalism [10]. 

Note that the probability of a member of the classical ensemble not satisfying the Lorentz 
gauge condition d^A^ = is zero by assumption, i.e., the Lorentz gauge is equivalent to 
the condition that the product (<9 M A M )P[A[ vanishes for all physical fields. For the quantum 
ensemble to satisfy this condition, as per the guiding principle discussed at the end of Section 
III.B above, one equivalently requires, noting Eq. (38), that 

(d^MA] = 0. 

As is well known, this constraint, if initially satisfied, is satisfied for all times [17] (as is the 
alternative weaker constraint that only the 4-divergence of the positive frequency part of 
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the field vanishes [10]). Hence the evolution of the quantum ensemble is consistent with the 
Lorentz gauge. It would be of interest to derive the consistency of this constraint directly 
from the equations of motion, Eqs. (CI) and (C2), for P and S. 

Finally, it is well known that one can also obtain the classical equations of motion via an 
alternative Hamiltonian, by exploiting the degree of freedom left by the Lorentz gauge to 
remove a dynamical coordinate (corresponding to the longitudinal polarisation). In partic- 
ular, since d^A^ is invariant under — > + d^x f° r an y function x satisfying d v d v x = 0, 
one may completely fix the gauge in a given Lorentz frame by choosing x such that A = 0. 
One thus obtains, writing A^ = (A , A), the radiation gauge A = 0, V • A = 0. The 
classical equations of motion for A (d u d u A = 0), follow, for example, from the Hamiltonian 

H R [A, E] = (1/2) J dx (E • E + | V x A| 2 ), (41) 

where E denotes the conjugate momentum density. Here Hr corresponds to the standard 
Lagrangian L = —(1/4) / dxF^F^, with A = 0. This Hamiltonian once again has the 
quadratic form of Eq. (32), and hence the exact uncertainty approach yields the correspond- 
ing Schrodinger equation 

ih(d*/dt) = H R [A, -ih(5/5A)]V (42) 

for a quantum ensemble of electromagnetic fields in the radiation gauge (this is also the form 
of the Schrodinger equation obtained via the Schwinger-Tomonaga formalism [18]). Note 
that the momentum density E is in fact the electric field, and hence, in this case, the exact 
uncertainty approach corresponds to adding nonclassical fluctuations to the electric field 
components, with the fluctuation uncertainty determined by the uncertainty in the vector 
potential A. 

D. Example: Gravitational field 

The gravitational field is described, up to arbitrary coordinate transformations, by the 
metric tensor g = (g^). The corresponding invariant length may be decomposed as [19] 

ds 2 = g^dx^dx" = -(a 2 - {3 ■ f3)dt 2 + 2(3 i dx' l dt + %jdx l dx j , 

in terms of the lapse function a, the shift function (3, and the spatial 3-metric 7 = (7^). 
The equations of motion are the Einstein field equations, which follow from the Hamiltonian 
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functional [19] 

if [7, tt, a, /3] = /" rfxa^ G [7, tt] - 2 / dx (3^^ , (43) 

where 7r = (tt 1 - 7 ) denotes the momentum density conjugate to 7, \j denotes the covariant 
3-derivative, and the Hamiltonian density Tic is given by 

H G = (l/2)G ij - fc ,[7]7r i % w - 2 (%[ 7 ](det 7 ) 1/2 . (44) 

Here is the curvature scalar corresponding to 7^, and 

£^[7] = (7ifc7ji + luljk ~ 7ij7fc/)(det7) _1/2 . 

The Hamiltonian functional H corresponds to the standard Lagrangian L = 
J dx (— det gy/ 2 R[g], where the momenta ir° and tt 1 conjugate to a and $ respectively 
vanish identically [19]. However, the lack of dependence of H on 7r° and tt 1 is consistently 
maintained only if the rates of change of these momenta also vanish, i.e., noting Eq. (Bl) 
of Appendix B, only if the constraints [19] 

5H/5a = H G = 0, 5H/5(3i = -2vr ij ^ = (45) 

are satisfied. Thus the dynamics of the field is independent of a and (3, so that these functions 
may be fixed arbitrarily [19]. Moreover, these constraints immediately yield H = in Eq. 
(43), and hence the system is static, with no explicit time dependence [19]. 

It follows that in the Hamilton- Jacobi formulation of the equations of motion (see Ap- 
pendix B), the momentum potential S is independent of a, (3 and t. Noting that tt = SS/Sj 
in this formulation, Eqs. (45) therefore yield the corresponding constraints 

5S 5S dS f5S\ 

for S. As shown by Peres [20], a given functional ^[7] of the 3-metric is invariant under 
spatial coordinate transformations if and only if (SF/Sji^y = 0, and hence the fourth 
constraint in Eq. (46) is equivalent to the invariance of S under such transformations. This 
fourth constraint moreover implies that the second term in Eq. (43) may be dropped from 
the Hamiltonian, yielding the reduced Hamiltonian 

H G [y,7T, a] = J dxaHoil^] (47) 
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in the Hamiltonian-Jacobi formulation [20, 21]. 

For an ensemble of classical gravitational fields, the independence of the dynamics with 
respect to a, f3 and t implies that members of the ensemble are distinguishable only by 
their corresponding 3-metric 7. Moreover, it is natural to impose the additional geometric 
requirement that the ensemble is invariant under spatial coordinate transformations. One 
therefore has the constraints 

5P 5P dP ( 5P\ 

for the corresponding probability density functional Ppy], analogous to Eq. (46). The first 
two constraints imply that ensemble averages only involve integration over 7. 

Noting Eq. (44), the Hamiltonian H G in Eq. (47) has the quadratic form of Eq. 
(32). Hence the exact uncertainty approach is applicable, and immediately leads to the 
Schrodinger equation 

ihdV/dt = J dxaH G [l,-ih(5/5l)]V (49) 

for a quantum ensemble of gravitational fields, as per the Corollary of Sec. III.B. 

As discussed at the end of Sec. III.B, we follow the guiding principle that all constraints 
imposed on the classical ensemble should be carried over to corresponding constraints on 
the quantum ensemble. Thus, from Eqs. (46) and (48) we require that P and S, and 
hence \l/ in Eq. (38), are independent of a, (3 and t and invariant under spatial coordinate 
transformations, i.e., 

{6V\ 

fa=w = 8r =0, Mr (50) 

Applying the first and third of these constraints to Eq. (49) immediately yields, via Eq. (44), 
the reduced Schrodinger equation 

Hah, -ih(6/6>y)]* = (-tf /2)^G m [l\-^l - 2 ( 3 )i?[ 7 ](det 7 ) 1/2 ^ = 0, (51) 

which may be recognised as the Wheeler-deWitt equation in the metric representation [19]. 
The consistency of Eqs. (50) and (51) is well known [19]. 

A notable feature of Eq. (51) is that the Wheeler-deWitt equation has not only been 
derived from an exact uncertainty principle: it has, as a consequence of Eq. (37), been 
derived with a precisely defined operator ordering (with Gijki sandwiched between the two 
functional derivatives). Thus the exact uncertainty approach does not admit any ambiguity 
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in the description of quantum gravity, unlike the standard approach [19]. Such removal of 
ambiguity is essential to making definite physical predictions, and hence may be regarded 
as an advantage of the exact uncertainty approach. 

For example, Kontoleon and Wiltshire [22] have pointed out that Vilenkin's prediction of 
inflation in minisuperspace, from a corresponding Wheeler- de Witt equation with "tunneling" 
boundary conditions [5], depends critically upon the operator ordering used. In particular, 
considering the class of orderings defined by an integer power p, with corresponding Wheeler- 
deWitt equation [5] 

r Ffl n f> 1 f> 2 

^ = (52) 



d 2 p d Id 2 TT/ 
+ "n 9 Trr-? U(a, 0) 



da 2 a da a 2 d(f) 2 

(for a Friedmann- Robertson- Walker metric coupled to a scalar field 0), Kontoleon and Wilt- 
shire show that Vilenkin's approach fails for orderings with p > 1 [22]. Moreover, they 
suggest that the only natural ordering is in fact the "Laplacian" ordering corresponding to 
p — 1, which has been justified on geometric grounds by Hawking and Page [23]. 

However, noting that the relevant Hamiltonian functional in Eq. (2.7) of Ref. [5] is 
quadratic in the momentum densities of the metric and the scalar field, the exact uncertainty 
approach may be applied, and yields the Wheeler-deWitt equation corresponding to p = —1 
in Eq. (52). Hence the criticism in Ref. [22] is avoided. One also has the nice feature that 
the associated Wheeler-deWitt equation can be exactly solved for this "exact uncertainty" 
ordering [5]. 



IV. DISCUSSION 



The two main results of this paper are (i) an exact uncertainty relation, Eq. (21), valid for 
all bosonic quantum fields, and (ii) the derivation of the quantum equation of motion, Eq. 
(37), from an exact uncertainty principle, for fields with Hamiltonian functionals quadratic 
in the momentum density. These results are conceptually connected by the notion of the 
momentum density splitting into the sum of classical and nonclassical components, with 
the uncertainty of the nonclassical component determined by the uncertainty of the field. 
Moreover, they generalise similar results obtained in Refs. [1, 2, 3, 4] for quantum particles. 

It is important to emphasise that the exact uncertainty approach in Sec. III.B does not 
assume the existence of a complex amplitude functional ^f[f], nor the representation of 
fields by operators, nor the existence of a universal constant H with units of action, nor 
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the existence of a linear operator equation in some Hilbert space. Only the assumptions 
of causality, independence, invariance and exact uncertainty are required, all formulated 
in terms of a single nonclassical element (the uncertainty introduced by the momentum 
fluctuation N). Since uncertainty is at the conceptual core of quantum mechanics, this is 
an elegant and pleasing result. 

The assumptions used also provide an intuitive picture for the origin of the Schrodinger 
equation for bosonic fields, as arising from nonclassical fluctuations of the momentum den- 
sity. Of course this picture has limitations - the fluctuations essentially arise from the 
uncertainty of the field itself, rather than from some external source, and hence are most 
certainly "nonclassical" rather than "classical" in nature. 

A minimalist interpretation of the exact uncertainty approach, based on Eqs. (33) and 
(36), is that every physical field has an intrinsic uncertainty, the nature of which precludes a 
deterministic relationship between the field and its conjugate momentum density. However, 
the degree of indeterminism in this relationship is precisely quantifiable, in a statistical sense, 
and is directly connected to the ensemble representation of the inherent field uncertainty. 

The exact uncertainty approach is very different to the "causal interpretation" of Bohm 
and co-workers [6] . In the latter it is assumed that there is a pre-existing complex amplitude 
functional ^[f] = y/P exp(iS/h) obeying a Schrodinger equation, which acts upon a single 
classical field via the addition of a "quantum potential", Q[P], to the classical Hamiltonian. 
It is further assumed that the momentum density is precisely g = 5S/5f, and that physical 
ensembles of fields have probability density functional P = \^>\ 2 . In contrast, the exact 
uncertainty approach does not postulate the existence of adjunct amplitudes and potentials, 
the Schrodinger equation directly represents the evolution of an ensemble rather than of 
an external amplitude functional acting on individual systems (and is derived rather than 
postulated), and the basic tenet in Eq. (33) is that g ^ SS/8f. 

It is of interest to consider the scope and limitations of the exact uncertainty approach to 
physical systems. Its applicability to nonrelativistic quantum particles [3, 4, 24], and now 
to bosonic quantum fields with Hamiltonians quadratic in the momentum density (which 
include all relativistic integer-spin fields), has been demonstrated. It is also, indirectly, 
applicable to the non-quadratic Hamiltonian functional of a nonrelativistic boson field (cor- 
responding to second quantisation of the usual nonrelativistic Schrodinger equation), in the 
sense that this field may be obtained as a low-energy limit of the complex Klein-Gordon 
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field [8] (to which the exact uncertainty approach directly applies). 

In this paper the basic Schrodinger equation for bosonic fields has been obtained, with 
the advantageous feature of a unique operator ordering in cases where other approaches are 
ambiguous. It would be of interest to consider further issues, such as the representation 
of general physical observables by operators (addressed for the case of particles in Ref. 
[3]), boundary conditions, infinities, etc, from the new perspective on the conceptual and 
technical basis of quantisation offered by the exact uncertainty approach. 

Finally, a major question to be addressed in the future is whether the exact uncertainty 
approach is applicable to fermionic quantum fields. These have two features which present 
challenges: the corresponding ensemble Hamiltonian is usually linear in the momentum 
density, and the anticommutation relations make it difficult to connect the equations of 
motion with corresponding classical equations of motion in the limit as H — > 0. One possible 
approach is to determine whether exact uncertainty relations exist for such fields, analogous 
to Eq. (21), as these would presumably suggest the statistical properties required by suitable 
"nonclassical" fluctuations. 

APPENDIX A: FUNCTIONAL DERIVATIVES AND INTEGRALS 

Here necessary definitions and properties of functionals are noted, including variational 
properties of functional integrals. Since manipulations in the paper are of a purely formal 
character, we avoid discussions of regularisations and discretizations needed to fully define 
certain quantities. 

A functional, F[f], is a mapping from a set of physical fields (assumed to form a vector 
space) to the real or complex numbers, and the functional derivative of F[f] is defined via 
the variation of F with respect to /, i.e., 



for arbitrary infinitesimal variations Sf. Thus the functional derivative is a field density, 
SF/Sf, having the value 8F/8f x at position x. For curved spaces one may explicitly include 
a volume element in the integral, thus redefining the functional derivative by a multiplicative 
function of x; however, this is merely a matter of taste and will not be adopted here. The 
functional derivative is assumed to always exist for the functionals in this paper. 




(Al) 
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It follows directly from Eq. (Al) that the functional derivative satisfies product and 
chain rules analogous to ordinary differentiation. The choice F[f] = f x > in Eq. (Al) yields 
Sf x i/Sf x = 5(x — x') as required in Eqs. (3) and (4). Moreover, if the field depends on some 
parameter, t say, then choosing 5f x = f x (t + St) — f x (t) in Eq. (Al) yields 

dt ~ at + J sf x at { ] 

for the rate of change of F with respect to t. 

Functional integrals correspond to integration of functionals over the vector space of 
physical fields (or equivalence classes thereof). The only property we require for this paper 
is the existence of a measure Df on this vector space which is translation invariant, i.e., 
/ Df = / Df for any translation f — f + h (which follows immediately, for example, from 
the discretisation approach to functional integration [8]). In particular, this property implies 
the useful result 

r fiF t 
j D f— = for J DfF[f] <oo, (A3) 

which is used repeatedly below and in the text. Eq. (A3) follows by noting that the finiteness 
condition and translation invariance imply 

= f Df (F[f + $f] ~ F[f]) = J dxSf x (| DfSF/Sf x 

for arbitrary infinitesimal translations. 

Thus, for example, if F[f] has a finite expectation value with respect to some probability 
density functional P[/], then Eq. (A3) yields the "integration by parts" formula 



J DfP(5F/5f) = -J Df(5P/5f)F 



Moreover, from Eq. (A3) the total probability, / Df P, is conserved for any probability flow 
satisfying a continuity equation of the form 

?L + J ix A_ lPVi] = , (A4) 

providing that the average flow rate, (V x ), is finite. 
Finally, consider a functional integral of the form 

I[F]= jDft(F,8F/8f), (A5) 
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where £ denotes any function of some functional F and its functional derivative. Variation 
of I[F] with respect to F then gives, to first order, 



51 = I[F + 5F] - I[F] = J Df {(d£/dF)6F + J dx [d^/d(5F/5f x )] [6(6F)/6f x ]} 

= jDf j(ae/0F) - / dxJj-[d£/d{6F/6f x )]} 5F 
+ Jdx J Df-^-{[d£/d(5F/5f x )]5F}. 



5 

Assuming that the functional integral of the expression in curly brackets in the last term is 
finite, this term vanishes from Eq. (A3), yielding the result 
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* (A6) 



analogous to Eq. (Al), where the variational derivative 51 /5F is defined by 

^L = K_ \dx — 

5F OF J 8 f x 

A similar result holds for multicomponent fields, with summation over the discrete index a 
in the second term. 



d(5F/5f x )\ 



APPENDIX B: HAMILTON- J ACOBI FIELD THEORY 

The salient aspects of the Hamilton- Jacobi formulation of classical field theory [9] are 
collected here, with particular attention to the origin of the associated continuity equation 
for ensembles of classical fields, required in Sec. III. 

Two classical fields /, g are canonically conjugate if there is a Hamiltonian functional 
H[f,g, t] such that 

df/dt = 5H/5g, dg/dt = -5H/5f. (Bl) 

These equations follow from the action principle 5A = 0, with action functional A = 
f dt [-H + / dx g x (df x /dt)}. The rate of change of an arbitrary functional G[f,g,t] follows 
from Eqs. (A2) and (Bl) as 

dG _ dG r (5G5H _ 5G 5H\ _ dG 
dt dt J \5f x 5g x 5g x 5f x J dt 

where { , } is a generalised Poisson bracket. 

26 



A canonical transformation maps /, g and H to /', g' and H', such that the equations of 
motion for the latter retain the canonical form of Eq. (Bl). Equating the variations of the 
corresponding actions A and A' to zero, it follows that all physical trajectories must satisfy 

-H + J dx g x (df x /dt) = -H' + J dx g' x (df' x /dt) + dF/dt 

for some "generating functional" F. Now, any two of the fields f,g,f',g' determine the 
remaining two fields for a given canonical transformation. Choosing / and g' as the two 
independent fields, defining the new generating functional G[f,g',t] = F + J dx f x g' x , and 
using Eq. (A2), then yields 



H> = H + f + fdx 



dt \6f x » X J dt \5g' x 

for all physical trajectories. The terms in round brackets therefore vanish identically, yielding 
the generating relations 

H' = H + dG/dt, g = 5G/5f, f' = 6G/Sg'. (B2) 

A canonical transformation is thus completely specified by the associated generating func- 
tional G. 

To obtain the Hamilton- Jacobi formulation of the equations of motion, consider a canon- 
ical transformation to fields /', g' which are time-independent (eg, to the fields / and g 
at some fixed time to)- From Eq. (Bl) one may choose the corresponding Hamiltonian 
H' = without loss of generality, and hence from Eq. (B2) the momentum density and the 
associated generating functional S are specified by the functional equations 

9 = jj, ^ + H[f,5S/5f,t]=0. (B3) 

The latter is the desired Hamilton- Jacobi equation. Solving this equation for S is equivalent 
to solving Eqs. (Bl) for / and g. 

Note that along a physical trajectory one has g' = constant, and hence from Eqs. (A2) 
and (B3) that 

dS dS r , 5S df x TT f , df x dA 
—r = t; — t- / dx -zir-^- = —H + dx q x —r- = — — • 
dt dt J 5f x dt J y dt dt 

Thus the Hamilton- Jacobi functional S is equal to the action functional A, up to an additive 

constant. This relation underlies the connection between the derivation of the Hamilton- 
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Jacobi equation from a particular type of canonical transformation, as above, and the deriva- 
tion from a particular type of variation of the action, as per the Schwinger-Tomonaga formal- 
ism [18, 25]. In the latter case the time parameter t is replaced by the multi-time parameter 
a. 

The Hamilton- Jacobi formulation has the interesting feature that once S is specified, 
the momentum density is determined by the relation g = 5S/5f, i.e., it is a functional of /. 
Thus, unlike the Hamiltonian formulation of Eqs. (Bl), an ensemble of fields is specified by a 
probability density functional P[f], not by a phase space density functional P[f,g}. In both 
cases, the equation of motion for P corresponds to the conservation of probability, i.e., to a 
continuity equation as per Eq. (A4). In particular, since in the Hamilton- Jacobi formulation 
the rate of change of the field / follows from Eqs. (Bl) and (B3) as the functional 



the associated continuity equation for an ensemble of fields follows from Eq. (A4) as [26] 



Eqs. (B3) and (B4) generalise immediately to multicomponent fields, and form the basis of 
the classical starting point in the derivation of the quantum equations of motion for bosonic 
fields in Sec. III. 

APPENDIX C: PROOFS OF THE THEOREM AND COROLLARY IN SEC. Ill 

Proof of Theorem (Eq. 36): From the causality and exact uncertainty assumptions 
one has Cov x (iV) = a(P, 8P/5f x , f x ). To avoid issues of regularisation, it is convenient to 
consider a position-dependent canonical transformation, f x — > A~ x f x , such that A [A] := 
exp[J dx In | det A x |] is finite. Then the probability density functional P and the measure 
Df transform as P — > AP and Df — > A~ l Df respectively, and so the invariance assumption 
requires that 



where u a and w a denote the vectors 5P/5f® and respectively, for a given value of x. 
Since A x can remain the same at a given point x while varying elsewhere, this homogeneity 



V x [f] = df x /dt = (5H/Sg x ) g=ss /sf , 




(B4) 



a(AP, AA^u, A^w) = A^a(P, u, w)A. 
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condition must hold for A and A x independently. Thus, choosing to be the identity 
matrix at some point x, one has a(AP, Au,w) = a(P,u,w) for all A, implying that a can 
involve P only via the combination v := u/P. The homogeneity condition for a therefore 
reduces to 

a(A T v, A^w) = A T a(v,w)A . 

Note that this equation is linear, and invariant under multiplication of a by any function of 
the scalar J := v T w. Moreover, it may easily be checked that if o and r are solutions, then 
so arc ot~ 1 o and tu^ 1 t. Choosing the two independent solutions a = vv T , r = (ww T )~ 1 , 
it follows that the general solution has the form 

a(v, w) = (3(J)vv T + 'j(J)(ww T )~ 1 

for arbitrary functions (5 and 7. For P = P1P2 one finds v = (vi,v 2 ), w = (wi,w 2 ), 
where the subscripts label corresponding subsystem quantities, and hence the independence 
assumption reduces to the requirements 

P(Ji + J2)Vivl + 7 (Jx + J 2 )(w 1 wJ)- 1 = A( JOuiuf + 7i(^i)(wiO _1 , 

P(Jl + J2)V 2 V2 + 7(Jl + J^iw^lY 1 = P 2 (J 2 )V 2 V2 + l2{J2){w 2 wlY l , 

for the respective subsystem covariance matrices. Thus j3 = Pi = f3 2 = C, 7 = 71 = 72 = D 
for universal (i.e., system-independent) constants C and D, yielding the general form 

[Cov x (N)] ah = C(6P/6f a x )(5P/5f b x )/P 2 + DW: h [f] 

for the fluctuation covariance matrix, where W x [/] denotes the inverse of the matrix with 
a6-coeflicient f x f x - Note that the latter term is purely a functional of /, and hence merely 
contributes a classical additive potential term to the ensemble Hamiltonian of Eq. (34). It 
thus has no nonclassical role, and can be absorbed directly into the classical potential (V) 
(indeed, for fields with more than one component this term is singularly ill-defined, and 
hence can be discarded on physical grounds). Thus we may take D = without loss of 
generality. Finally, the positivity of C follows from the positivity of the covariance matrix 
Cov x (N), and the theorem is proved. 

Proof of Corollary (Eq. 37): First, the equations of motion corresponding to the 
ensemble Hamiltonian H' follow via the theorem and Eqs. (31) as: (a) the continuity equation 
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Eq. (29) as before (since the additional term does not depend on S), which from Eq. (32) 
has the explicit form 




and (b) the modified Hamilton- Jacobi equation 

dS/dt = -SH'/SP = -H[f, 5S/5f, t] - 5(H' - H)/5P. 
Calculating the last term via Eq. (36) and Eq. (A6) of Appendix A, this simplifies to 

t+WSItM-K^J* (^ + ^) =0. (02) 

Second, writing ^ = P 1 / 2 exp(iS/h), multiplying each side of Eq. (37) on the left by ^P -1 , 
and expanding, gives a complex equation for P and S. The imaginary part is just the 
continuity equation of Eq. (CI), and the real part is the modified Hamilton- Jacobi equation 
of Eq. (C2) above, providing that one identifies C with h 2 /A. 
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